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RESUMEN 


El aluminio es uno de los metales de mayor uso industrial y domestico debido a su baja densidad, alta conductividad 
térmica y eléctrica, es liviano y tiene buenas propiedades mecánicas. Este elemento se obtiene de la alúmina, siendo 
la temperatura de fusión de la alúmina 2040%C; esta es demasiada elevada para poder abordar el proceso 
industrialmente. Por lo que, aplicando energía eléctrica por kilogramo de alúmina [6.3 -14] Kwh/KgAl, a través de un 
método de electrólisis en una celda de criolita (NazAlF6ó) fundida, se baja la temperatura del proceso a 1000%C. La 
gran corriente eléctrica y la alta temperatura en la celda traen como consecuencia, entre otros, desgate de ánodos, 
formación de costras en las paredes de la celda, efecto anódico no deseado, contaminación e inestabilidad en el 
control del proceso. Con el fin de disminuir el problema, abordamos el método numérico de diferencias finitas para 
resolver la ecuación de transferencia de calor asociadas al sistema físico de la celda. Para esto, estudiamos los perfiles 
de temperaturas en el espacio y el tiempo mediante las herramientas computacionales en el simulador PDETOOL de 
MATLAB y el manipulador y simulador matemático MAPLE. Finalmente, los resultados analíticos y numéricos 
alcanzan una buena aproximación convergente. 


PALABRAS CLAVE: Reducción de aluminio, Transferencia de calor, Diferencia finita, Pruebas de convergencia y 
estabilidad. 


ABSTRACT 


Aluminum is one of the largest industrial metals and domestic use because of its low density, high thermal and 
electrical conductivity, is lightweight and has good mechanical properties. This element is obtained from the alumina, 
being the melting temperature of alumina 2040*C; This is too high to approach the process industrially. So, applying 
electrical energy per kg of alumina [6.3 -14] Kwh / kgal, through a method of electrolysis in a cell of cryolite 
(Na3AIF6) melt, the process temperature is lowered to 1000 * C. The high power and high temperature in the cell 
brought as a consequence, among others, desgate anode, crusting on the walls of the cell, anode effect unwanted 
pollution and instability in the process control. To lessen the problem, we address the finite numerical method for 
solving the heat transfer equation related to physical differences in the cell system. For this study the temperature 
profiles in space and time using computational tools in MATLAB PDETOOL simulator and the manipulator and 
MAPLE mathematical simulator. Finally, analytical and numerical results reach a good convergent approach. 


Keywords: Reduction of aluminum, heat transfer, finite difference, convergence and stability tests. 
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1. INTRODUCCIÓN 

El aluminio se produce a través de un proceso de reducción electrolítica en una celda. Los 
componentes principales que reaccionan en la celda son el polvo de alúmina u óxido de 
aluminio (Al203), junto a un catalizador formado por un baño químico donde su principal 
componente es la criolita (Na3A1F6), en menor cantidad se agregan otros compuestos para 
mejorar el proceso (Peter, 2013). A través de los ánodos se hace pasar una corriente de 
320KAmp, la tensión anódica aplicada con respecto al cátodo es de 4.5V, esto transforma la 
alúmina sólida (polvo de alúmina) en aluminio líquido a un rango de temperatura promedio de 
960—C. La técnica más usada consiste en el proceso de reducción elaborado por los 
investigadores Paul Héroult y Charles Hall en 1886. En la figura 1, se muestra este proceso de 
reducción electrolítica y sus principales componentes (Peter, 2013). 


Carbon anodes 


Figura 1. Diagrama del proceso Hall-Héroult 


La principal reacción química que ocurre en este proceso se puede representar por la 
siguiente ecuación 


241,0 5h IC sólido) > AAL eiguido) els 3CO Gas) (1) 


3 (Disuelta) 


Las soluciones que se han planteado para modelar los problemas presentados en esta 
reacción electroquímica (1) tienen un alto grado de complejidad debido a la cantidad de 
variables no independientes que intervienen en este proceso. Las variables son la corriente 
eléctrica, la temperatura, campo electromagnético, la presión de gases, entre otros. Todo el 
proceso de transferencia de calor está caracterizado por la ecuación no lineal de difusión de 
calor en el espacio (Gutierez y Troyani, 1999), (Gutierez y Troyani, 2004). 


LT y pco (2) 
al ó) 9Ll y) al ót 


Donde T es la temperatura, q es la generación de calor, k la conductividad térmica, C calor 
específico, t el tiempo y pes la densidad del material. 


Este artículo presenta en la parte I y II la descripción del proceso de reducción de aluminio, 
los modelos y las ecuaciones analíticas que lo rigen y se discretizan las ecuaciones en 1D. La 
sección Ill, describe el esquema del seudocódigo para 1D en lenguaje Maple. La parte IV y V se 
generan las simulaciones en Pdetool y Maple. Finalmente, en VI se analizan los resultados y en 
VII las conclusiones. 
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HN. ECUACIONES CARACTERÍSTICAS DEL PROCESO 


El aluminio en su fase líquida dentro de la celda se puede aproximar a una cuba tipo ladrillo a 
fin de aplicar métodos numéricos y obtener una solución aproximada [2]. La figura 2 es un corte 
transversal de esta celda donde se muestran la dimensión y sus principales elementos: ánodo, 
cátodo, baño electrolítico (criolita) y el aluminio líquido. 


Barra 


metálica 
4 TT conductora 


(Electrodo +) 


Baño Electrolítico ll 


Figura 2. Corte transversal de una celda de reducción 


Las ecuaciones (1) y (2) dirigen del proceso químico y térmico respectivamente que se lleva 
a cabo en la celda. No obstante, la fórmula (2) de difusión necesita el parámetro voltaje- 
corriente de generación de calor. 


f= 0. .[VV 1 (3) 


Donde o. es la conductividad eléctrica y V el voltaje entre el ánodo y cátodo de la celda 
electrolítica. 


2.1 Esquema de diferencia finita en 2D 


Para el desarrollo de este trabajo consideramos a la conductividad térmica k de la ec. (2) 
constante. Siendo esto así, se puede escribir de la siguiente forma 


k (e) g _ OT 


peplóx? Oy” pep Ot 


OT TY q _oT 
7 LEA (4) 
Ox?  0y pes 0% 
donde a =k/pCp. 
Los términos en segunda derivada, 0T/0x? y 0*T/0y?, de la ec. (4) se pueden desarrollar 
tomando aproximaciones en la serie de Taylor [4], con dependencia en una sola coordenada, de 
la forma 


TIT Tra0o rana, 
3 ] (5) 
7) +01) 
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T(x,)=7G1)+T (10M + ray E + 
(6) 


TU) Sl +0(M%). 


Aquí (), () y (") son la Ira, 2da y 3ra derivada respectivamente. Si truncamos la serie hasta 
segundo orden, obtenemos: 


_ Tx.) -2T(x)+T(%,_,) 


E +0(45). (0) 


T(x,) 


El término O(h?) es el error de truncamiento de la serie. Ahora bien, incorporamos la ec.(7) 
como el esquema de discretización de la ec.(4), la cual permite escribir las segundas derivadas 
aproximadas (o de diferencia finitas), en el punto nodal ¡, ;, de la forma 


o"T Toja 27, + Ta 2 
2|— ET(Ax”) (8) 

2 2 

ox” |, Ax 

y 

ST] (Ta -2T,+Ta, 

7 ATA), (9) 
Y lo y 


Donde 4x y 4y definen las dimensiones de cada nodo y los subíndices i y j sirven para designar 
las posiciones x y y de los puntos nodales discretos. La aproximación en diferencia finitas para 
la derivada respecto el tiempo en la ec.(4) se expresa de la forma 


OT 
Ot 


n+l n 
2 Ta — UN 
At 


(10) 


Dj 


Donde 41 es el intervalo de tiempo entre cálculos sucesivos, el superíndice n se utiliza para 
denotar la dependencia con respecto al tiempo de 7, y la derivada con respecto al tiempo se 
expresa en términos de la diferencia en temperatura asociada con los tiempos nuevo (n+1) y 
anterior (n). Esto se muestra gráficamente en la figura 3. 


i-L,¡e«= 


i,j-1 i+L j 


Figura 3. Modelo cúbico de celda para representar los nodos en los tiempos n y (n+1) 


Haciendo O = ek y sustituyendo las ecuaciones. (8), (9) y (10) en (4), se obtiene la forma 
explícita de la ecuación en diferencia finitas para el nodo interior i, j 
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ae 53 ZE, + LA ] Os; 
Ay PCp 


n+ n oLAt n n n 
Tp ha PN A N E 
a (12) 
A At 
Ed 541 li =21,, +1, sl La, 
Ay PCp 


Esta última ecuación es explícita, pues las temperaturas nodales desconocidas para el tiempo 
nuevo se determinan de manera exclusiva mediante temperaturas nodales conocidas en el 
tiempo anterior. Es decir, la ecuación (12) proporciona la distribución transitoria de temperatura 
al avanzar en el tiempo con el uso de intervalos de 4f. 


En el caso del esquema de diferencias finitas implícito, se puede obtener una reducción en la 
cantidad de tiempo. Para esto, se evalúan todas las otras temperaturas en el nuevo tiempo (n+1), 
en lugar del tiempo anterior (n). Se considera entonces que la ecuación (4) proporciona una 
aproximación en diferencias hacia atrás para la derivada con respecto al tiempo. Entonces, 
reordenando nuevamente en (12), distribuimos los términos n+1 en el lado izquierdo y los 
términos con n en el lado derecho 


a LEI E e pS 


y 14, 
n 
Ar (13) 
n+l n+l rn O; 
SL a SL ja = E +A —É 
PCp 


Donde por comodidad se ha tomado s,=a.4t/4x? y sy=0a4t/Ay?. De la ec.(13) es evidente que la 
temperatura nueva del nodo ¡,j depende de las temperaturas nuevas de sus nodos contiguos que, 
en general, se desconocen. Por tanto, para determinar las temperaturas nodales desconocidas en 
t + At, las correspondientes ecuaciones nodales deben resolverse simultáneamente. 


2.2 Definición de convergencia y estabilidad 
Se dice que una solución numérica converge a la solución analítica cuando para cada punto, en 
el espacio de las variables independientes, la primera tiende a la segunda al refinar la grilla de 
cálculo, es decir 
e (14) 


At>0 
Ax>0 


x; fijo 

t” fijo 
Nótese que al efectuar el paso al límite indicado en la Ec. (14) los índices ¡ y n deben variar 
(específicamente; j,n —> 0) para que x; y 1” permanezcan fijos (por ejemplo, x; =jAx y t" = 
nAt). 


La estabilidad de la solución numérica respecto de la analítica significa acotar, de alguna 
manera, la diferencia [T(x;, t”) = TF"|. En problemas de valores iniciales, que resultan en 
algoritmos iterativos, puede adoptarse la siguiente definición [4] que se establece si, 


< 00 as) 
n>0 
Ax,AL fijo 


ENT 
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es decir, la diferencia entre la solución analítica y numérica para cada nodo del eje espacial, 
permanece acotada cuando se avanza indefinidamente el cálculo en la dirección temporal; 
manteniendo fijos los intervalos de discretización. Una restricción de este tipo sólo puede 
proveer, en principio, condiciones necesarias para la estabilidad. 


2.3 Solución analítica por el método de separación de variables. 

En esta sección desarrollamos un método sencillo para buscar una solución analítica en 2D de 
la ecuación (4). En este caso, consideramos q = 0 y en condiciones de isotropía espacial para 
un plano de lados £ y H. Es decir en coordenadas cartesianas la ec.(4) se transforma en 


2 2 
El E O<x<L y 0<y<H (16) 


+ E > 

ox oy” Ot 
Entonces, proponemos el método de separación de variables [4], de la forma 
T(x, y,t) =v(x, y)u(1) (17) 


Sustituyendo la ec.(16) en (17), obtenemos 


2 2 
> 4 (18) 


not vlóx oy 
donde A es la constante de separación. Luego se obtienen de las dos ecuaciones 
p+Au=0 (9) 
a v+4v=0 (20) 
aquí (') significa la derivada con respecto a £ y el operador V es la derivada con respecto a las 


coordenadas espaciales x y y. Luego, por condiciones iniciales y de fronteras de Neumann y 
Dirichlet, escribimos 


1(0)= 4, V(0,y)=V,0, V(L, y) =V,, ,v (5,0) = Vio y v(x, H) = Vo: 
Luego de (19), obtenemos la solución en el tiempo 


u=e”, (21) 


y para la parte espacial aplicamos nuevamente el método de separación de variables, 
v(x, y) = X (Y (y). (22) 


Sustituyendo la ec.(22) en (20), obtenemos 
—+24=-— =P (23) 


donde £ es la constante de separación. Luego, de.(23), obtenemos 


AE: co 42) + Casen 4) (Q4) 
o o 
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Y =C, cos 2) +C, sol 2) (5) 
10 10 


donde Cy, C», C3, C4 son constantes y Ó es la constante de separación, la cual tiene la relación 
9 =4-— f. Sustituyendo las expresiones (23), (24) en (22) y conjuntamente la expresión (21) 
en (17), se llega a la solución general: 


TIE uer[Ceof 22) + a) 
a a 


(26) 
(Geo) +C se +) 
o o 


Ahora bien, si la ecuación (26) se desea evaluar en una dimensión, hacemos d¿=0 y y=0, para 


obtener 
T(x,t) = [ÉS + E] (27) 
a a 


Una solución más simple se obtiene al hacer C2=0 


a 


Donde el producto 7,C> representa una sola constante con unidades de temperatura. La 
expresión (28) es la solución analítica para ser implementada en la simulación de la ec. (2) en 
1D. 


T(x,t) = ES (28) 


2.4. Solución analítica en 1D utilizando el manipulador MAPLE. 


Otra solución analítica unidimensional se ha obtenido utilizando el lenguaje de programación 
matemático MAPLE. 


Enter a PDE. 


0 o a 
or T(x,t) dá TG0) 0 


> Solve PDE 
(T(x,t) = _ FIG) _F2(0) £ Where 
E Fl)=_c, Flo, ES F2AY=a c, ra) 
> PDEtools [build] 
a it ASE 
Ta) =e*_ Ce a O a (29) 


donde C,, C> y C3 son constantes. 


TI. ESQUEMA DEL SEUDOCÓDIGO 


A continuación se presenta un seudocódigo para resolver la ecuación de transferencia de calor 
en 1D, con O =0 considerando el método de diferencias finitas, conjuntamente con la 
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asignación del dato inicial y las condiciones de fronteras mixtas de Neumamn y Dirichlet. Para 
esto mostramos detalles de la construcción de un seudocódigo en al manipulador y simulador 
matemático MAPLE. 

- Seudocódigo en MAPLE para la Ecuación de Transferencia de Calor. 


+ Datos fisicos y químicos de la ecuación de calor (en la reducción de aluminio) 


k=237: + Conductividad térmica (W/m.K) 
p=2702: + Densidad del material (Kg/m"3) 
Cp» =903: + Calor específico (J/Kg.K) 

k 
194 + Parámetro de la ecuación de calor 


a 


++ Entrada de datos que define el tiempo 


f, =0: ++Tiempo inicial 
tp =1: Tiempo final 
At=0.001: Tiempos sucesivos de cálculos 
Í == 
ntiempo = dl Tiempo en el espacio 


++ Entrada de datos que definen el espacio 
x =0: *+Punto inicial del alambre 


x=: 
dk 
Ax=WVAx-d” : HEspaciado computacional 


FPunto final del alambre 


Xy Xp ] , 
; + Número de puntos espaciales 


nespacio = 


++ Condiciones de frontera: Neumann-Dirichlet 
for j from 0 to ntiempo do 


HOME Le? 
T[L, /]:= Te *0os(E) 


od: 
+ El dato inicial 
for i from 0 to nespacio do 


T[1,0] := T,cos(2%%) : 
10 
od: 


4 El programa 
for j from 1 to ntiempo-1 do 
for i from 1 to nespacio-1 do 
T[1,¡+1]:=7[i, ¡]+ 


At-a 
2 
x 


(T[¡+1 7]-2-7[2, ¡]+7[-1, 1D) 


od. 
od : 
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Iv. MODELAJE Y SIMULACIÓN EN PDETOOL 


La simulación en PDETOOL de MATLAB se realizó con los parámetros del aluminio dados 
en (Cengel, 2004; Incropera, 1996) las tablas 1, 2, 3 y 4 de características internas de los 
materiales el INCROPERA, usando un corte transversal de la celda de reducción tipo rectángulo 
laminar 2D de aluminio. 


Tabla 1. Parámetros de la ec. de transferencia de calor. 


Siglas | Descripción Magnitud | Unidades 
rho | Densidad 2072 Kg/m? 

C Capacidad calórica 903 JKgK 
K Coef. de transf. de calor 237 w/meK 

Q Origen de calor 0 “C 

h Coef. de conv. de transf. calor | 1 

Text _ | Temp. Externa 40 “C 


Tabla 2. Condiciones de borde superior (Dirichlet). 


Siglas | Descripción Magnitud | Unidades 
g Flujo de calor 0 

q Coef. de transf. de calor 0 w/m*K 

h Coef. de conv. de transf. calor | 1 

r Aporte de calor 960 “C 

Tabla 3. Condiciones de borde inferior (Dirichlet). 
Siglas | Descripción Magnitud | Unidades 
g Flujo de calor 0 
q Coef. de transf. de calor 0 w/mK 
h Coef. de conv. de transf. calor | 1 
r Aporte de calor 660 “C 
Tabla 4. Condiciones de borde laterales (Neumann). 

Siglas | Descripción Magnitud Unidades 

g Flujo de calor 0 

q Coef. de transf. de calor | 0 w/m*K 

h Coef. de convecc. de 1 

transf. de calor 
r Aporte de calor 0 “C 


En las figuras 3 y 4 se muestra la variación de temperatura entre 10 y 50 segundos. En las 
figuras 5 y 6, luego de transcurrir 500 segundos, se observa la uniformidad en los rangos de 
temperatura de 950*C hasta 700*C, el aluminio se encuentra en su fase liquida. 


Figura 3. Simulación PDETOOL a los 10 segundos. 
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6 


Figura 5. Simulación PDETOOL a los 500 segundos. 


Figura 6. Simulación PDETOOL a los 1000 segundos. 
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n £ 


Figura 7. Perfil de temperatura en periodo permanente. 


Nodo 41 
Nodo 64 
Nodo 89 
Nodo 69 


Nodo 41 
Nodo 64 
Nodo 83 
Nodo 69 


050 10D 150 200 250 300 350 400 450 500 0100 200 300 400 500 600 700 800 900 1000 


Figura 8. Variación de temperatura en nodos centrales. 


Las curvas de comportamiento térmico en los nodos centrales de la placa de aluminio líquido 
indican que la temperatura se estabiliza. En la figura 8, se observa que la estabilidad en los 
nodos 41 y 69 se obtiene después de los 400 segundos de reacción en la celda electrolítica. 


Todas estas gráficas son el resultado de la simulación del módulo de tratamiento de 
ecuaciones diferenciales PDETOOL de Matlab y muestran el comportamiento de una sección 
transversal o lámina rectangular de la celda de reducción de aluminio Hall-Héroult. 


V. MODELAJE Y SIMULACIÓN EN 1D CON MAPLE 


En esta sección partimos del hecho de que la solución numérica difiere de la solución analítica 
debido a dos causas: el error de truncamiento y los errores de redondeo. Además, se asume que 
la solución numérica es convergente al afinar la malla de cálculo y disminuir el error de 
truncamiento, mientras los errores de redondeo permanecen bajo control. En el estudio 
numérico y analítico de la ecuación de transferencia de calor en 1D se implementó un algoritmo 
numérico en MAPLE, bajo las condiciones iniciales y de frontera de Dirichlet y Neumann: 


- Frontera con la solución de la ecuación (26) 


TQ.) =T,e", TL =T,e*cos(2L) y 


T(c0)=7, .cos(É%) (30) 


- Frontera con la solución (27) 
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TOD=CCE "OO" y 


TN =C Cer, ae 
> 123 E CL 
T(x%,0) =C,C,e 4 * + 3 (31) 
e 


Con el fin de validar los resultados numéricos y analíticos, seleccionamos Ax = 0.018 y un 
At = 0.001; los cuales mejoraron el refinamiento de los puntos en el tiempo. 


En las figuras 9, 10 y 11 se muestra la solución numérica con línea de puntos negros y la 
solución analítica en línea sólida de la ec. (31), del mismo modo, las figuras 12, 13 y 14 detallan 
el valor absoluto de la diferencia entre las dos soluciones planteadas. En estas gráficas se 
muestra la convergencia de la solución numérica con la solución analítica, la cual indica que el 
error global de truncamiento es muy pequeño para los tiempos de (0, 10 y 100 segundos 
respectivamente). 


intercepción de la Solución Numérica y Analítica en t=0(5) 
700 


Figura 9.Solución analítica £ numérica a 0 seg. 


Intercepción de la Solución Numérica y Analítica ent=10(5) 
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Figura 10. Solución analítica 8 numérica a 10 seg. 


Ciencia e Ingeniería 2015 Vol. (2) N* (1) A. Oliver, E. Lobo, H. Martínez y R Acuña 
EVALUACIÓN MEDIANTE DIFERENCIA FINITA DE TRANSFERENCIA DE CALOR 
QUE OCURREN EN EL PROCESO DE REDUCCIÓN DE ALUMINIO 


Intercepción de la Solución Numérica y Analífica en t=100(5) 
700 


Figura 11. Solución analítica £ numérica a 100 seg. 


Diferencia |[T(mm-Tfanal)| en t=0(5) 
1,x10:7 


3,x10:8 


6,x10"* 


[Téxm)Tfanal)| 


4x108 


Figura 12. Diferencia de analítica y numérica a los 0 seg. 


Diferencia |Tímem-Tíanal | en t=10(5) 
40 


[Téumm)Tfanal )| 20 


Figura 13. Diferencia de analítica y numérica a los 10 seg. 
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Diferencia |Tímum-Tfanal )| en t=100(5) 


404 


[Tun Tfanal )| 20, 


Figura 14. Diferencia de analítica y numérica a los 100 seg 


Con el fin de ampliar sobre los errores que describen la convergencia y estabilidad de las 
soluciones de la ecuación de transferencia de calor, se prueban las condiciones iniciales y de 
frontera de la ec. (30). 


La implementación del método de separación de variables se describe gráficamente en las 
figuras 15, 16 y 17, donde se muestra la solución analítica de la ecuación de transferencia de 
calor con línea sólida y la solución numérica en puntos negros para una corrida de 0, 10 y 100 
segundos respectivamente. Del mismo modo, en las gráficas 18, 19 y 20 se muestra el valor 
absoluto de la diferencia entre la solución analítica y la solución numérica para los tiempos de 0, 
10 y 100 segundos respectivamente Aquí se puede ver una gran aproximación de la solución 
numérica con la analítica. 


Intercepción de la Solución Numérica y Analífica en t=0(5) 
300 


Figura 15. Método de separación de variables 0 seg. 
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Infercepción de la Solución Numérica y Analítica en t=10(5) 


Figura 16. Método de separación de variables 10 seg. 


Intercepción de la Solución Numérica y Analífica en t=100(5) 


Figura 17. Método de separación de variables 100 seg. 


Diferencia |T(men-Tíanal )| en t=0(5) 
1,x10:? 


3,x10:8 


6,x10"* 


[Té Tfanal )| 


4x108 


Figura 18. Separación de variables |An — num] 0 seg. 


Ciencia e Ingeniería 2015 Vol. (2) N* (1) A. Oliver, E. Lobo, H. Martínez y R Acuña 
EVALUACIÓN MEDIANTE DIFERENCIA FINITA DE TRANSFERENCIA DE CALOR 
QUE OCURREN EN EL PROCESO DE REDUCCIÓN DE ALUMINIO 


Diferencia |T(mum-Tíanal)| en t=10(5) 


[Tun Tfanal )| 
D,02 


Figura 19. Separación de variables |4n — num] 10 seg 


Diferencia |T(nmum-Tíanal )| en t=100(5) 
] 


[Téuem)Tfanal)| 


Figura 20. Separación de variables |4n — num] 100 seg. 


El problema de estabilidad de un algoritmo consiste en analizar la sensibilidad de los 
resultados a pequeñas variaciones en los datos de entrada. En términos laxos, se denomina 
estable a aquél cuyos resultados muestran poca sensibilidad a esas variaciones y viceversa. 
Obviamente, estas consideraciones requieren ser expresadas con mayor rigor matemático. Con 
Fundamento en lo expuesto anteriormente, se observa que, para ambos casos, la diferencia 
|An — num] en t=0seg, es muy pequeña (107), este valor puede ser tomado como un cero. 


Para hacer esta evaluación se somete a prueba tomando la ecuación (31) y variando el 
parámetro de la constante C/. Primero se evalúa con 0.15 en espacio tiempo y luego, se reduce 
la constate a 0.015, es decir, diez veces menos. Las gráficas de la figura 21 muestran que la 
solución numérica para 7(x,t) cambia drásticamente ante las variaciones de C.. 


Ciencia e Ingeniería 2015 Vol. (2) N* (1) A. Oliver, E. Lobo, H. Martínez y R Acuña 
EVALUACIÓN MEDIANTE DIFERENCIA FINITA DE TRANSFERENCIA DE CALOR 
QUE OCURREN EN EL PROCESO DE REDUCCIÓN DE ALUMINIO 


C,=0.15seg C,=0.015seg 
Figura 21. Temperatura en la coordenada espacial y el tiempo para C:=0.15 y 0.015 


VI. ANÁLISIS DE LOS RESULTADOS 


En las gráficas de las figuras 9, 10, 11, 12, 13 y 14; se muestra la convergencia de la solución 
numérica con la solución analítica de la ecuación (31), la cual indica que el error global de 
truncamiento es muy pequeño. Por lo tanto, se puede confirmar que la solución numérica 
converge a la solución analítica. Este resultado muestra una línea de material que se enfría en un 
corto tiempo y hacia el centro de la línea y luego, comienza a calentarse. No obstante, las 
figuras 15, 16, 17, 18, 19 y 20 muestran una buena convergencia, pero, el material tiende a 
enfriarse en el tiempo. 


Hay dos posibilidades para determinar la causa de la inestabilidad, presente en la figura 21. La 
primera es que sea una característica del problema que puede reflejar, por ejemplo, una 
inestabilidad física, en cuyo caso se trata de un problema mal condicionado o matemáticamente 
inestable; en nuestro caso, la comparación del comportamiento de la temperatura para el valor 
de la constante C;=0.15, con el comportamiento, para C¡;=0.015, denota una característica del 
problema. La segunda posibilidad, más común que la anterior, es que el algoritmo esté 
pobremente construido; entonces, se estaría en el caso de un algoritmo mal condicionado o 
numéricamente inestable. 


VII. CONCLUSIONES 


Estimamos que para disminuir el error de truncamiento, se recalcule los resultados con el 
truncado efectuado en un nivel superior (por ejemplo, incluyendo más términos en la evaluación 
de una serie). En este caso desarrollar fórmulas que den cuenta, a menos, de la dependencia de 
los errores de truncamiento respecto de los parámetros del problema, por ejemplo, el número de 
términos utilizado en la evaluación de una serie. 


La simulación en PDETOOL facilita la visualización del intercambio de calor en las 
diferentes capas que componen el proceso de reducción de aluminio, gráficamente se observa la 
interacción entre los contornos que aportan o no aportan calor al proceso de reducción de 
aluminio, en los periodos transitorios y permanentes. 


La ecuación diferencial utilizada para simular el proceso de reducción de aluminio, nos 
muestra que es convergente en el periodo permanente, por mostrar la temperatura en forma 
constante con un diferencial pequeño en el tiempo. Se puede concluir que la validación de las 
soluciones numéricas y analíticas depende de los datos asignados a las constantes de integración 
y al modelo desarrollado. 
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